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Abstract 

We present a fully quantum-electrodynamical formalism suitable to evaluate the spontaneous 
emission rate and pattern from a dipole embedded in a non-absorbing and lossless multilayer 
dielectric structure. In the model here developed the electromagnetic field is quantized by a proper 
choice of a complete and orthonormal set of classical spatial modes, which consists of guided 
and radiative (partially and fully) states. In particular, by choosing a set of radiative states 
characterized by a single outgoing component, we get rid of the problem related to the quantum 
interference between different outgoing modes, which arises when the standard radiative basis is 
used to calculate spontaneous emission patterns. After the derivation of the local density of states, 
the analytical expressions for the emission rates are obtained within the framework of perturbation 
theory. First we apply our model to realistic Silicon-based structures such as a single Silicon/ air 
interface and a Silicon waveguide in both the symmetric and asymmetric configurations. Then, 
we focus on the analysis of the spontaneous emission process in a silicon-on-insulator (SOI) Slot 
waveguide (a 6 layers model structure) doped with Er'^"'" ions (emitting at the telecom wavelength). 
In this latter case we find a very good agreement with the experimental evidence [M. Galli et ai, 
Appl. Phys. Lett. 89, 241114 (2006)] of an enhanced TM/TE photoluminescence signal. Hence, 
this model is relevant to study the spontaneous emission in Silicon-based multilayer structures 
which nowadays play a fundamental role for the development of highly integrated multifunctional 
devices. 

PACS numbers: 78.67.Pt, 78.20.Bh, 42.50.Ct, 42.70.Qs 
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I. INTRODUCTION 



It is well known that the environment surrounding an excited atom affects its rate of 
spontaneous emission (SE): enhanced SE in a resonant cavity was first predicted in the 
pioneering work by Purcell [1] and, later on, an inhibited SE in a small cavity was shown 
by Kleppner |2]. Such an effect can be explained either by classical electromagnetism, in 
terms of a self-driven dipole due to the reflected field at the dipole position, or in the 
framework of quantum electrodynamics, as emission stimulated by zero-point fluctuations 
of the electromagnetic field. As long as the coupling between the atom and the field is weak, 
both descriptions yield the same results [5]. In such a weak coupling regime, the SE rate 
can be calculated within first order perturbation theory by applying the Fermi's Golden 
Rule, and is proportional to the local coupling of the atomic dipole moment to the allowed 
photon modes, i.e. to the local density of states (LDoS) [U [5]. The modifications in the 
electromagnetic boundary conditions induced by the surrounding material alter the density 
of states as well as the SE rate: when the LDoS vanishes, then the SE process is inhibited, 
while, when an increase in the density of states occurs, the rate of SE can be enhanced over 
the free space value. 

A large amount of work, both theoretical and experimental, has been devoted to the 
analysis of the SE from emitters (such as atoms, molecules or electron-hole pairs) embedded 
in dielectric environments of varying complexity. In an homogeneous medium with dielectric 
constant e, it has been shown by Glauber [B] that the SE rate relative to the free space 
value, is enhanced when e > 1 and reduced for £ < 1, as it has been demonstrated also 
experimentally by Yablonovitch [7]. In spite of, or rather, thanks to its relative simplicity, 
the single interface has been subject of a constant research [H [9l [101 El [121 [131 El [ISl [161 [IZ] 
which still goes on, since it is the ideal system where experimental and theoretical analysis 
can be performed in order to get the basic understanding necessary to investigate more 
complex structures. The double interface has also been widely studied, especially as a 
waveguide slab, i.e., an high-index core surrounded by low- index cladding layers, with both 
a quantum electrodynamical [IHl [19l [201 1211 [22] and classical [23j approach. In systems 
characterized by more than two interfaces, both the technology and the theory needed 
become more demanding, but the expected effects turn out to be also more interesting. For 
instance, among multiple dielectric layers structures, planar microcavities have been subject 
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of intense research in last years [211 [251 ISSl |27], due to their abihty to considerably affect 
the density of states and thus strongly modify the emission into a particular mode, which is 
of crucial importance for the development of light emitting devices. 

In this work we study the SE rate in a non-dispersive and lossless multilayer dielectric 
structure by applying a fully quantum electrodynamical formalism. With respect to previous 
published works, which generally deal with a specific dielectric structure, our main aim 
is to develop a model suitable for more than one configuration, thus taking into account 
all the possible modes (and the related SE rates) which can be excited in the examined 
structure. While our discussion tackles the problem of the spontaneous emission from a 
theoretical point of view, the results derived can provide a useful quantitative insight into 
the modifications of the atomic radiative processes which occur in realistic structures. As an 
example, we apply our method to evaluate the SE rate in silicon-based optical waveguides, 
which nowadays can be tailored according to different geometries, from simple waveguides 
(three-layers geometry) to multilayers-like configurations. These structures are characterized 
by a high index contrast and are able to confine and guide the light very efficiently in 
nanometer-size spatial regions as a result of total internal reflection. The waveguiding and 
confining properties, together with the low propagation losses (typical of Silicon) and the 
good compatibility with complementary metal oxide semiconductor technology, make them 
very attractive for the future development of highly integrated multifunctional optoelectronic 
and photonic devices (see e.g., Refs. [2HI EHl EO])- Furthermore, with an embedded optically 
active layer, these structures can also be exploited as silicon-based optical sources. With this 
purpose, one of the most promising configurations recently proposed is the Slot waveguide 
[3T| [32] : this system consists of a thin (few nanometers thick) layer (slot) of low-index 
material bounded by two high-index material regions (typically Silicon), which are the core of 
an optical waveguide; the high-index contrast interfaces at the slot are able to concentrate the 
electromagnetic field in very narrow spatial regions (<^ A), thus leading to an enhancement in 
the radiation-matter interaction. A theoretical investigation [33J of the emission properties 
of a slot waveguide doped with Erbium ions and embedded in optoelectronic devices, as well 
as the experimental evidence [31] of an enhanced light-matter coupling, have been already 
presented, but a full quantum mechanical analysis of the spontaneous emission processes in 
this kind of structures is still lacking. Here, we face this problem by applying the developed 
formalism to evaluate the SE rate of a dipole embedded first in a single Silicon/air interface. 
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then in both a symmetric (high-index contrast) and an asymmetric Sihcon waveguide, and 
finally in a Slot waveguide. 

In order to build up a quantum electrodynamical theory of the SE process, the elec- 
tromagnetic field must be first decomposed into the normal modes supported by dielectric 
structure under consideration. This is needed in order to set up a second quantized form 
of the electromagnetic field, and then to express the local density of states and SE rate by 
application of Fermi's Golden Rule. The LDoS can be also derived within a quantum electro- 
dynamic and Green's function formalism as often done in the literature, using either a scalar 
or a dyadic Green's function, see e.g. Refs. [101 HU EH ES]- In a generic multilayer structure, 
the normal set of modes, i.e., a complete and orthonormal set of solutions of Maxwell equa- 
tions for the considered structure, is well known [37] and consists of a continuous spectrum 
of radiative modes and a discrete one composed of guided modes, defined for both transverse 
electric (TE) and transverse magnetic (TM) polarizations [38]. Guided modes are trapped 
by the highest refractive index layer (if any), and are evanescent in both half spaces - the 
lower and upper cladding - surrounding the multilayer structure. Radiative modes can be 
either fully or partially radiative. The former, similar to free space modes, extend over the 
whole space and propagate out of the dielectric structure from both cladding layers as out- 
going plane waves, while partially radiative modes propagate from the cladding layer with 
higher refractive index only, being evanescent (due to total internal reflection) along the 
lower refractive index cladding. The modes, found as the elementary solutions of Maxwell 
equations with proper boundary conditions, have more than one representation, since one 
needs to characterize the asymptotic behavior of the radiative states, such a characterization 
being not unique. The standard set of radiative modes, originally introduced by Carniglia 
and Mandel [39], and which is generally applied to describe the interaction of a radiating 
system with the electromagnetic field in a dielectric structure, is not very convenient for SE 
analysis though. In this paper, we chose to apply a set of radiative modes characterized by 
a single outgoing component only. Such a choice leads to a simple definition of the LDoS for 
radiative states, avoiding the difficulties related to the treatment of the interference between 
different outgoing modes (see Refs. [iQlHT]), which arise when the standard set of radiative 
modes based on the triplet incident-refiected-transmitted waves is used. Furthermore, the 
emission rates in the lower and upper half-spaces of a generic multilayer structures - or, in 
general, the SE patterns - can be easily calculated. 
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The paper is organized as follows. In Sec. II the field modes supported by multilayer 
dielectric structure are listed and described. We show that the basis of radiative states which 
has been used for the quantization of the electromagnetic field in the considered dielectric 
structure, can be obtained from the standard set of radiative modes by a Time-Reversal 
transformation. In Sec. Ill we perform a standard quantization of the electromagnetic field, 
and in Sec. IV a second quantized form for the atom-field interaction term of the whole 
system Hamiltonian is set up and then used (in the electric dipole approximation) to derive 
the expressions of the LDoS and the SE rate as a function of the dipole position. In Sec. V 
the spatial dependence of the SE rate will be examined for several structures of interest. A 
short summary of the results is given in Sec. VI. 

II. SYSTEM AND FIELD MODES 

The system we are investigating is depicted in Fig. [T] it is made up of M dielectric layers 
(stack) which are parallel to the xy plane and assumed to be infinite along the x and y 
directions. Each layer is dj (j = 1,...,M) thick and the surrounding media, i.e., the lower 
(layer 0) and the upper (layer M+1) claddings, are taken to be semi-infinite. Each of the 
M + 2 media is supposed to be lossless, isotropic, and homogeneous along the vertical {z) 
direction. Hence, the dielectric constant £:(r) = e(p, z) is a piecewise constant function in 
the z direction and it will be denoted as Sj = ej{z) in each of the M + 2 dielectric media. 



In order to develop a quantum theory for the spontaneous emission of a dipole embedded 
in such a dielectric structure, the classical electromagnetic modes, which are needed in the 



expansion of the electromagnetic field operators (see Sec. Ill), must be first specified. The 
modes are found as the solutions of the following eigenvalue problem 
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which results from the homogenous Maxwell equations for the electric and magnetic fields 
E, H having harmonic time dependence exp(— icut), and with the condition V ■ H = being 
fulfilled. The set of these fundamental modes is complete and orthonormal. 



H(r) 



(2) 




FIG. 1: Schematic view of the multilayer dielectric structure. The lower (layer 0) and the upper 
(layer M + 1) claddings with dielectric constants Eq and eM+i) respectively, are taken to be semi- 
infinite and surround the stack made by M dielectric layers, each one having a thickness of dj and 
characterized by an average dielectric constant ej, j = 1, M. 

the orthonormality condition being expressed by 

J H;(r)-H,(r)dr = 5^,. (3) 

c 

The electric field eigenmodes, which can be obtained from E(r) = i — ^ H(r), are also 

uje{r) 

orthonormal according to the following condition [39] : 

J £(r)E;(r)-E,(r)dr = 5^,. (4) 

Since the whole system is homogenous in the xy plane the field modes will be factorized 
as E(r,t) [H(r,t)] = e'''^'+'^rPE{z) [Il{z)], where ky = A;||A;|| = {k^,ky) is the in-plane 
propagation vector. 

In a lossless multilayer dielectric structure, the complete set of orthonormal modes consists 
of an infinite number of radiative modes and a finite number of guided modes. The former 
can be classified into two types. Fully radiative modes, akin to free-space modes, radiate 
in both the lower and upper cladding, while partially radiative modes radiate only in the 
cladding with the higher refractive index, propagating out of the smaller index cladding as 
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evanescent waves with exponentially decreasing amplitude. Guided modes propagate along 
the dielectric planes only, being trapped (confined) by the highest refractive index layer and 
characterized by an evanescent field profile in both claddings. 

Whereas the guided modes are completely specified by the Maxwell equations and the 
proper continuity conditions across the dielectric boundaries, the radiative modes are not, 
and their asymptotic behavior at infinity (when z ±00) has to be characterized. Such a 
characterization, however, is not unique. 




FIG. 2: The radiative modes [38] in a multilayer dielectric structure, (a): the standard set of 
modes based on the triplets incident-reflected-transmitted waves, {ii, ri,ti} for waves incoming 
from the lower cladding, {iu, ru,tu} for waves incoming from the upper cladding, (b): the set 
of modes specified by a single outgoing component and two incoming (towards the stack) waves, 
{Xq , Wo, Xm+i} for states outgoing in the lower cladding, {Wm+i , Wq, Xm+i} for states outgoing 
in the upper cladding. The notation refers to TE-polarized modes; for TM polarization one needs 
the replacements W ^ Y and X ^ Z. 



The standard choice for radiative modes [see Fig. [2|^a)] assumes one incident wave incoming 
towards the stack of M layers either from the lower or from the upper cladding, and two 
outgoing waves, one being refiected (on the same side of the incoming one) from the stack 
and the other being transmitted (on the opposite side) across it. This set of modes, originally 
introduced by Carniglia and Mandel [39] for the quantization of the electromagnetic field in 
a dielectric interface, is orthonormal and complete [12] and it has been widely employed to 
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characterize the radiative states in structures hke dielectric waveguides [T8l [21] and planar 
dielectric microcavities [23, ES] ■ 

Such a choice, however, is not the most convenient when dealing with radiation emission 
analysis. As shown in Fig.[2]^a), both the reflected and the transmitted components (the pairs 
{ri, tu} and {ru, ti}), which belong to two different modes, contribute to the total emission 
in a given direction. As pointed out by Zakowicz |10], the computation of the radiative 
density of states turns out to be problematic since the quantum interference between the 
two different outgoing modes has to be explicitly taken into account. Thus, interference 
terms must be considered when emission in either the upper or the lower layer is evaluated. 
As shown in the reply by Glauber and Lewenstein [41j, interference effects cancel out only 
when the emission in the upper cladding layer ({ru, ti}) and in the lower one ({ri, tu}) are 
combined. To avoid this subtle interference problem, and to be able to calculate the fraction 
of emission in either the upper or the lower cladding layer, a more suitable way to define the 
radiative states in emission problems is to choose a set of modes based on a single outgoing 
wave. This component comes together with two incoming waves propagating towards the 
structure, as shown in Fig. |2](b) with the triplets {Wm+i , Wq, Xm+i} and {Xq , Wq, Xm+i} 
for states outgoing in the upper and lower cladding layers, respectively. By using this mode 
decomposition, the total emission signal is thus completely specified by one outgoing mode 
only - either by the component labeled as Wm+i for radiative states outgoing in the upper 
cladding, or by the component Xq for states outgoing in the lower cladding - and interference 
terms never arise. 

It is worth to notice that this set of radiative modes can be obtained from the stan- 
dard one previously discussed (which is specified in terms of the incoming waves) after 
application of Time-Reversal (TR) operator. A Time-Reversal operation transforms modes 
propagating along the positive 2;-direction into modes propagating along the negative one, 
and therefore a mode incoming from the upper (lower) layer, into a mode outgoing from 
the upper (lower) layer. Hence, as from Figs. |2](a)-(b), the radiative modes {Xq , Wq, Xm+i} 
and {Wm+1 , Wq, Xm+i} are the TR counterparts of the triplets {ii, ri,ti} and {iu, ru,tu}, 
respectively. Furthermore, since the algebraic properties are invariant under Time-Reversal 
operations, also the new set of radiative states is orthonormal and complete. The following 
rule for the Time-Reversal operation over a generic spatial mode oc e*"^^ propagating along 
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the 2;-direction with wavevector g, can be estabhshed: 

^iqz TR^ ^-iq*z _ 

The transformation rule Eq. ([s]) apphes to both fully and partially radiative modes: outgoing 
modes towards positive (negative) z being oc e*''^ (oc e~^^'^), specified by the real wavevector 
g, keep the plane-wave-like character, turning into outgoing modes towards negative (pos- 
itive ) z oc e"*"^^ (oc e*''^). The evanescent modes oc e""^'^', characterized by the imaginary 
wavevector q = in, keep the exponentially decaying profile after the transformation Eq. (|5]). 

A mode decomposition characterized by a single outgoing component, like the one here 
described, has been already used to specify the radiative modes in dielectric interfaces [I3l [16] 
as well as in slab waveguides [19]. It has also been used in the formally similar problem of 
diffraction losses in photonic crystal waveguides |13]. Here, we extend its application within 
a quantum electrodynamical theory suitable to the analysis of radiation emission in generic 
multilayer structures. 

In the following a detailed description of both guided and radiative profiles is given. 

A. Radiative modes 

As previously introduced, the set of radiative states consists of a single outgoing com- 
ponent propagating outward from the whole structure and two other waves propagating 
towards it. In each of the M layers the field is a superposition of two modes propagating in 
opposite directions (with respect to the z-direction) . The modes are specified and labeled 
by the propagation wavevector k = (k||,g), where the 2;-component q in each of the M-l-2 
media, is given by 

q, = \fl^^l, i = 0,...,M + l. (6) 

Let us denote by = z ^ k\\ the unit vector which is orthogonal to both ky = A;||A;|| and 
z and set zi = —di/2, zj = zj-i + dj^i with j = 2, ...,M + 1. With the implicit time 
dependence e"*"^*, the field profiles for TE polarization are given by 

El^{p,z) = '^^i^^^E^^{k^^,z), (7) 



H^^(P,;.) = '^^-[Hl^{^,z)z + H{^{^,z)^], (8) 
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where V is a normalization box- volume which disappears in the final results. The expressions 
for the amplitudes E"^^, Hj^ and Hj"^ as well as the method used to obtain them are detailed 
in Appendix |X} 

For TM-polarized radiative modes the field profiles are given by: 

H^f(p,.) = ^.ek„i/™(A:||,.), (9) 

Er(P'^) = '^^ — [ET'{^,z)z + E{^i^,z)k], (10) 



where H'^^, E™ and E7^ are given in Appendix A. 



B. Guided modes 



In order for the whole dielectric structure to support a set of guided modes, (at least) one 
of the dielectric constant ej {j = 1, M) of the M inner layers has to fulfill the constraint 

£j = ^max > £0 , £M+1 , (H) 

The guided modes, which are in-plane propagating and evanescent along the z direction, 



are labeled by the in-plane wavevector ky = k\ik\\ and the mode index a [a > 1 if Eq. (11) 
holds] in a joint single index /i = (ky, a). By we denote the 2;-component of the guided 
mode wavevector, 



0^2 



qj, = \/e,J-kl J = 1,...,M, (12) 

where o;^ = uj^^^a is the frequency of the a-th guided mode. In the upper (j = M + 1) 
and lower (j = 0) claddings qj ^ is purely imaginary, qj^ = iXjfi where Xi/x ~ ~ ^ ^ 
and hence the mode field oc exp(±igj^z) decays exponentially along the z direction. In the 
following we give their explicit form, which results from a generalization of the standard 
waveguide field modes (see e.g. Refs. HI]). 

The guided modes for TE polarization are given by: 

E5f(p,.) = ^^^4„i?™(A:||,.), (13) 
alfiP,^) = ^[i/r(A:|,,^)i + ifP(A:||,^)^ll], (14) 
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where S* is a normalization surface which cancels in the final results, and E"^^, Hj^ and 



H 



TM, 



if™ are given in Appendix |B For TM polarization, the guided modes have the following 
field profiles: 

(15) 
(16) 



S 

TM 



where H'^^, E™ and E7^ are given in Appendix 



B 



III. FIELD QUANTIZATION 

In this Section the canonical quantization of the electromagnetic field in a non-uniform 
isotropic dielectric medium described by a piecewise constant permeability e{r) is per- 
formed [HI Hni US]- The electric displacement vector and the magnetic induction (a unit 
magnetic permeability is assumed) are simply given by the relations 

D = £(r)E, B = H. (17) 

The starting point is the quantization of the vector potential A which is defined by the 
familiar relations 

B = V X A , (18) 
1 f)A 

E = _v$--^. (19) 

c dt ^ ^ 

We use here the generalized Coulomb gauge [0] defined, in absence of external charges, by 
the choice $ = and the relation 

V-[£(r)A] = 0, (20) 

which automatically satisfies the transversality condition on D, V ■ D = V ■ [£:(r)A] =0, 
and is consistent with the equation of motion for the vector potential A 

Vx(VxA) + ^^=0. (21) 

In order to obtain a second-quantized Hamiltonian for the free photon field, we first introduce 
the classical Hamiltonian function Tiem, i-e. the total electromagnetic energy, 

Hem= ^^n(r,t)A(r,t)dr-£ 

= i^|^[.(r)E(r)2 + B(r)2]dr, (22) 

12 



where is a quantization volume, 11 = e(r)A(r, t)/c^ is the canonical momentum, and 
L = — / [£:(r)E(r)^ — B(r)]dr is the Lagrangian function from which Eq. (21 ) follows after 

Hamilton's principle. The vector field operator A is then expanded in normal modes: 



A = ^(27r;i^kn)'/'[aknAk„(r)e-^">^"* + aLAL(r)e 



(23) 



where a|^^ (flkn) are Bose creation (destruction) operators of field quanta with energies hio-^n 
satisfying the usual commutation relations 



[Skn, Sk/^/] — (5k,k'^n,n' 5 [Skn, ^k'n'] — [Sk„, a-k/„/] — , 



(24) 



n being a generic index labeling the corresponding eigenmode characterized by the wavevec- 
tor k. 



/ e(r)AL(r)-Ak.„,(r)dr 



2~^k,k'^n,n' , 
kn 



(25) 



From Eq. ( 25 ) the following orthonormality conditions [HI EHl SS] for the electric and mag- 
netic fields follow: 

(26) 



/ £(r)E*„(r) ■ Ek'„Kr)dr = 5k,k'5. 
Jv 

Bk„(r) • Bk'n'(l")dr = 5k,k'5n,n' 



n,n' ) 



(27) 



V 



Finally, from Eqs. (22) and (23), one gets the well known second-quantized form for the free 
photon field: 

" ' (28) 



k,n ^ ^ 



IV. EMISSION RATES 



In this Section the spontaneous transition rate of an excited atom embedded in a non- 
uniform dielectric medium is calculated. We suppose that the atom, located at position z 
and initially in the excited state |x) (having energy hu^) undergoes a spontaneous dipole 
transition to its ground state |g) (having energy hUg) thereby emitting a photon of energy 
huo = — The total Hamiltonian of such a system can be written as 



7-A ) 



(29) 



where is the free- field Hamiltonian given by Eq. (28), is the free-atom Hamiltonian, 



Ha = huJo\x){x\ + ;ia;g|g)(g| 



(30) 
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and H-y-A is the atom-field interaction term which, in the electric dipole approximation and 
near the atomic resonance a; ~ tuo, reads [17] as 

iJ^_A ~ (fT+d + a_d*) ■ E(r, t) , (31) 

where o"_ = |g)(x| and 0"+ = |x)(g| are the atomic down- and atomic up- transition operators, 
respectively, and d = dxg = (x|d|g) = |d|ed is the dipole matrix element, d = er being the 
atomic dipole operator of the atom located at r. The electric field operator E(r) can be 



obtained from the vector potential operator A through Eqs. (18) and (23), 

E(r,t)= i J](27r;k^kn)^/^[aknEk„(r)e 



k,n 

- aLEL(r)e^--*] • (32) 

We assume that the interaction between the excited two-level system and the electro- 
magnetic field in the dielectric medium is not too strong, so that the transition between two 
states can be studied within the framework of perturbation theory. Let us then consider the 
initial \i) and the final |/) states of the combined atom-radiation system: initially there are 
no photons and the atom is in the upper (excited) level, \i) = |0) |x); in the final state one 
photon is emitted in any mode of the electromagnetic field of frequency (-^kn and the atom 
is in the lower (ground) level, |/) = |lkn,) ® |g). According to Fermi's Golden Rule (see e.g. 
Ref. [38]) the spontaneous emission rate F = F(r) of an atom located at position r is 

m = %J2\^fA-A\^)f6{u,~Uf), (33) 
/ 

where hwi and huf are th energies of the initial and final state, respectively. By insertion 
of Eq. (32) in the above expression, and using the commutation rules for a^n and a|^„, the 



spontaneous decay rate finally reads 

4 ^ I d I ^ 

r(r) = |Ekn(r) ■ edP^kn5(^o - ^kn) • (34) 

k,n 

By taking into account the i-th cartesian component of the eigenmode Ek„(r), the 
contribution F, to the total emission rate can be written as 

ri(r) = 1^1 ^° J,(lc;o, r) , (35) 
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where Ji{u!o, r) is the i-th contribution to the local density of states (LDOS) [U |5] J^ujq, r): 

Ji^o^r) = ^Ji{uJo,r). (36) 

i 

In a multilayer dielectric structure, an excited dipole can decay either as a radiative or 



a guided eigenmode. As discussed in Sec. II A, the radiative modes are specified by the 
propagation vector (kj|,g) of the outgoing component. Hence, in Eq. (34), k = (k||,g) and 



n={p ,j) is a double index specifying the final state parameters, namely the field polarization 
p = TE, TM and the cladding layer j in which the emission occurs, j = for emission in 
the lower cladding and j = M + 1 for emission in the upper cladding. For what concerns 
the guided modes, k = ky and n={p,a), where a is the guided mode index introduced in 
Sec. II B Furthermore, since the dielectric function e{r) = e{z) = Sj is homogenous in each 



layer, the spontaneous emission rate will be expressed as a function of the z coordinate only. 

For both decay channels (radiative and guided) two contributions to the total emission 
rate can be distinguished: (i) the emission rate Fy due to the decay of horizontal dipoles, i.e. 
in-plane oriented dipoles {i^ = x or id = y), which couple to both TE- and TM-polarized 
fields, (ii) the rate F_l due to the decay of vertical dipoles (^d = ^) which interact with 

TM-polarized modes only. For randomly oriented dipoles, the total averaged emission rate 

2 1 

can thus be written as F = -Fii H — Fi . In the rest of this Section we derive the exact 

3 " 3 

expressions for the emission rates into both radiative and guided modes. 



A. Emission rates into radiative modes 



For each propagation wavevector k = (k||,g) the frequency = (fcjj + g^)^^^ of the 

radiative modes has to satisfy the relation 

^y<^l=^J>' J=0, M + l, (37) 
where ej = Eq (^m+i) if the emission occurs in the lower (upper) cladding. With tUkn = ^-y 



in Eq. (34), the emission rate into the radiative modes T = T(z) is thus given by 

rw = ^ E E EEK,(p.'-)--.i^ 

p=TE,TM j=0,M+l k|| q 

X uj^S{uJo — co^) , (38) 
15 



where the TE- and TM-polarized fields E™(p, 2;) and E™{p,z) are given by Eq. Q and 



(10), respectively. It is convenient to re-write the emission rate as a function of the LDoS 



for radiative states Jra.d{^o,z) according to 
whith 



^{z) = ^^1^1 ^° J^^j(u;o, z) , (39) 



p=TE, TM j=0, M+1 

X pj(k||,co')dk|| , (40) 

and pj(k||,co') being the one-dimensional (ID) photon DoS at a fixed in-plane wavevector ky, 
for radiative modes outgoing in the medium j: 

,2 ,,,2~ 



9 



^v^-o (41) 



27rc 









where L = V/S is the width of the normalization box in the ^-direction (which disappears 
in the final expression of the SE rate) and 6 [9(x) = 1 (= 0) if a; > (x < 0)] is the 



Heaviside function. It is worth to stress that, by using Eq. (40) with the basis of radiative 



states discussed in Sec. II A we get rid of any ambiguity in the definition of the LDoS: for 
each outgoing radiative mode (j = or j = M + 1) the LDoS is defined by a single mode- 
component only and thus any difficulty related to interference effects between components 



of different modes is avoided. Also, due to the Heaviside function in Eq. (41), emission into 
partially radiative modes occurs only in the cladding with the higher refractive index. 
From Eqs. (39)-(41), and after the introduction of spherical coordinates in the (k||,g) 



space, 

kj| = (fcjsin^cos0, kjsmesmcf)) , 0e[O, 2tt] , 0e[O,7r/2] , (42) 

the single contributions to the total emission rate T{z) due to the decay of horizontal and 
vertical dipoles are easily obtained: 

rrW = ^ E ef£^'\E^^{^^ = k,sine,z)\'sinede, (43) 



j=0,M+l 
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r™(.) 



^ nc[e(z)]2 , 



/■7r/2 

^ ef^ \E^^{k\\=kjsme,z)\^smede, (44) 

=0, M+l 

^ ' = kjsine, z)\hmede , (45) 

0,M+1 



where the field amphtudes E™, and Ej'^ are given by Eqs. nAlh, i\A6v and i\A5\\ 

respectively. 



B. Emission rates into guided modes 



According to Eq. (34), the spontaneous emission rate for the decay into guided modes 



having frequency Uk^^a is given by 

47r2|dP 



Viz) 



h 



(46) 



p=TE,TM a kii 



<^gui('^05 z) 



where the fields EZ^(p, 2;) and E?'^(p, 2;) are given by Eq. (13) and Eq. (16), respectively. 



the sum extends over all the a guided modes, and the 2D LDoS Jgui(i-o'o, z) is given by 

S 



(2vr) 



_rpT7i r-TAT „, J 



(47) 



p=TE,TM a 

The emission rates rjT^, TV^ and T™ can be easily obtained after integration over ky of 



Eq. (46): 



II 



rr(z) 



hujQ 

Idp27rc2 



2 ^ 



2 '^0 



^5:|E™(A:|, = A:o^.), 

^1 

a 
a 



71 TM/ 



(48) 
(49) 
(50) 



where E™, and Ej_^ are given by Eqs. nBlh, (^B7| and rtB6|, respectively. In the 



expressions given above, k^ = k^{u! = Uq) and Vq = ((ico'k||a/'^^||)wk||Q=i^o in-plane 
wavevector and the group velocity of the a— th guided mode calculated at the dipole emission 
frequency uo, respectively. The wavevectors k^ as functions of the frequencies can be found 
as the poles (which are real ones for guided modes) of the transmission amplitude t = I/T22 
of the whole dielectric structure, T being the total transfer matrix. 
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V. APPLICATIONS 



In this Section we apply the formahsm previously developed in order to investigate the SE 
process in realistic multilayer structures. As a typical high-index dielectric material, we take 
Silicon (nsi=3.48). After the analysis of a single Silicon/air interface, we will examine and 
compare the emission and confinement properties of different Silicon waveguides, namely 
a standard waveguide slab consisting of a Silicon core surrounded by two cladding layers 
with the same refractive index (symmetric configuration) or different ones (asymmetric con- 
figuration) and the silicon-on-insulator Slot waveguide. The SE rate has been evaluated 
for dipoles emitting at Aq = c^o/c = 1.55 /im which is the typical emission wavelength of 
Erbium ions (Er^+) often used as the active layer of Silicon-based light sources (see e.g. the 
review paper by Kenyon |19]). All the rates shown have been normalized with respect to 
the vacuum emission rate F = Fq = (4|dpu;Q)/(3/ic^) of a randomly oriented dipole. 



Figure |3] shows the normalized spontaneous emission F/Fq for a Silicon/air interface as func- 
tion of zjX. The emission rate for horizontal dipoles decaying into TE- and TM-polarized 
modes [see Fig. [sj^a)] varies continuously through the interface, as required from the continu- 
ity condition of the tangential field component at a dielectric boundary, while the emission 
rate for vertical dipoles (which couple to TM-polarized modes only), is discontinuous at 
the same point [see z = in Fig. [3|^b)] due to the discontinuity of the z component of the 
electric field. Far from the interface boundary, when zjX ^ 1, the spontaneous emission 
rate (for both horizontal and vertical dipoles) is scaled by the refractive index according to 
F(z) = e{z)Vis, in agreement with earlier works |6l |12l [13], with oscillations around the 
average value. The contributions to the total emission rate due to the decay into partially 
and fully radiative modes are shown in Figs. |3](c)-(d): in the dielectric half space the emis- 
sion is mainly due to the partially radiative states (see the thick solid lines for zjX < 0) 
which also characterize the profile of the total emission rate in the proximity of the interface 
boundary [see the thick solid lines at values z/X between and 0.2 in Figs. [3]^a)-(c) and in 
Figs. |3]^b)-(d)] and decay exponentially (in form of evanescent waves) in the free half space 
far from it. Hence, the evanescent component of partially radiative modes (which does not 
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FIG. 3: The normalized spontaneous emission rate as function of the dipole position zjX (A = Aq) 
in a dielectric/air interface. The half space z/A < is made by Silicon with refractive index 
nsi = 3.48. (a): the contribution from horizontal (in-plane oriented) dipoles decaying into TE- 
and TM-polarized modes, (b): the contribution from vertical {z oriented) dipoles which couple to 
TM-polarized modes only, (c) and (d): the total emission rates due to fully and partially radiative 
modes for horizontal and vertical dipoles. 
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FIG. 4: The spontaneous emission rates F = Fiow and F = Fup, for light outgoing in the lower and 
upper half-space, respectively, as function of the dipole position zjX for the same dielectric/air 
interface depicted in Fig. [3j The rates are normalized to the total emission rate Ftot = Tiow + Tup- 



contribute to the total energy flux and it is hidden in standard far field experiments) turns 
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out to be relevant in radiation emission analysis, since it strongly affects the radiative life- 
time r = 1/r in the vicinity of the interface boundary. Moreover, for vertical dipoles in a 
generic dielectric/air interface, one can analytically work out that, in the limit of a very large 
refractive index 1, the emission into partially radiative modes (which is the dominant 
one) at the discontinuous boundary, is given by T[z — > 0~)/ro = and T{z 0+) = n, 
in agreement with an earlier work by Loudon fl2\ [50] . 

The contributions to the total SE rate corresponding to light emitted either in the lower or 
in the upper layer are shown in Fig.|4]as a function of the dipole position. It is worth to notice 
that, within our model, these quantities are straightforwardly obtained by selecting the single 
outgoing radiative mode (see Fig. [2] in Sec. |ll]), through the index-layer j = (emission into 
the lower cladding) or j = M + 1 (emission into the upper cladding) in Eq. (40). These rates 
could also be obtained by using the standard basis with a single ingoing and two outgoing 
components [39], but in this case the interference terms between the two modes of Fig.2a 
must be explicitly calculated [10]. Thus, the present approach using the basis with a single 
outgoing component is especially useful for calculating radiative patterns and the emitted 
light in the lower/upper half spaces, which is a physically and technologically important 
problem for light emitting structures like LEDs and vertical laser diodes. 

The spontaneous emission rates for a symmetric Silicon waveguide are shown in Figs. [5] 
and|6]as functions of z/d, where d is the thickness of the Silicon core {\z\ < d/2) surrounded 
by air. Such a structure supports a finite number of guided modes and, since the upper 
{z > d/2) and lower {z < d/2) claddings have the same refractive index, only fully radiative 
modes can be excited and propagate out from the waveguide. By choosing a thickness 
d = Xq, one can calculate 7 TE and TM guided modes whose contribution to the total 
emission rate is significantly greater than the contribution due to the emission into radiative 
modes, as it can be seen by comparison of Figs. [5]^a)-(b) with Figs. |6](a)-(b). Furthermore, 
the emission rate in the core [see the spatial range —1/2 < z/d < 1/2 in Figs. [5]^a)-(b)] is 
close to the bulk value risiFo. 



The influence of an increasing number of guided modes is investigated in Fig. [7], where 
the emission rate has been evaluated as a function of the so called dimensionless photonic 
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FIG. 5: The normalized spontaneous emission rate into guided modes for a symmetric Silicon 
waveguide with air claddings as a function of z/d, d being the thickness of the Silicon core. The 
refractive index in the half spaces \z\ > d/2 is 1 and the core thickness has been taken equal to 
vacuum emission wavelength Aq. (a): the contribution of horizontal dipoles. (b): the contribution 
of vertical dipoles. 

thickness oj^d/c, while keeping the dipole position at the centre {z = 0) of the waveguide. 
It can be noticed that, with increasing thickness d, the contribution from the new arising 
modes is associated with the appearance of distinct features in the emission pattern such as 
dips and peaks. 



Furthermore, for vertical dipoles, the spontaneous emission rate is drastically suppressed for 
waveguide thickness below d = 0.5{c/ujq) = 0.5(Ao/27r) [see Fig. [7](a)], whereas, for the same 
range of thicknesses, the total emission rate from horizontal dipoles [see the thick solid line 
in Fig. ^h)] is F ^ 1.3 ^3.3 Fq and mainly due the excitation of TE guided modes. Also, for 
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FIG. 6: The normalized spontaneous emission rate into radiative modes for a symmetric Silicon 
waveguide with air claddings as a function of z/d, The refractive index in the half spaces \z\ > d/2 
is 1 and the core thickness has been taken equal to vacuum emission wavelength Aq. (a): the 
contribution of horizontal dipoles. (b): the contribution of vertical dipoles. 



thicknesses d > 2c/uJo = Aq/vt, the contributions to the total emission rate due to horizontal 
and vertical dipoles become comparable and close to the bulk value risiFo. 

The above results, which follow from the mode decomposition based on a single outgoing 



component for radiative states (see previous Sec. II A), are in agreement with those shown 



in earlier works [T8l EH, [22] and which have been obtained by using the standard set of 
radiative modes based on the triplets incident-reflected- transmitted waves. 



We now apply our model to study the SE in an asymmetric dielectric waveguide, i.e., a 
waveguide with different refractive indices in the lower and upper claddings (both values 
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FIG. 7: The normalized spontaneous emission rate for a symmetric Silicon waveguide with air 
claddings evaluated as a function of the photonic thickness (wq d) /c and for a dipole kept at the 
middle of the Silicon core, (a): the emission from horizontal dipoles. (b): the emission from vertical 
dipoles. 

being of course smaller than the core one). Due to the asymmetry, the condition for total 
internal reflection can be met for incidence angles beyond the limit one, and also partially 
radiative modes, which are evanescent in the lower index cladding, can thus be excited in 
such a structure. Figures |8] and [9] show the z dependence of the emission rates into guided 
and radiative modes, respectively, for an asymmetric waveguide made by a Silicon core 
bounded by a Silicon Oxide (Si02) lower cladding, and by air in the upper half space acting 
as upper cladding. For a thickness d = Xq, there are now 7 TE and 6 TM guided modes and 
the asymmetry-induced modiflcations in the emission pattern are clearly seen, especially 



for emission into radiative modes. Figure 10 shows the emission rate as a function of the 



core thickness: there are no guided modes for thicknesses smaller than d ^ 0A2c/u!q and 
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FIG. 8: The spatial dependence of the normahzed spontaneous emission rate into guided modes for 
the asymmetric Sihcon waveguide (Si02/Si/air). The thickness d of the SiUcon core (nsi = 3.48) 
is equal to the vacuum emission wavelength Aq. The lower cladding (z < d/2) is made by Si02 
("■3102 = 1.45) and the refractive index of the upper cladding {z > d/2) is 1. (a): the contribution 
of horizontal dipoles. (b): the contribution of vertical dipoles. 



the emission rate is thus sustained by radiative modes only (see the continuous thick line). 
For d > c/ujq = Ao/27r the emission is mainly due to guided modes, the contribution from 
TE polarized modes being larger. However, as a consequence of slab asymmetry leading 
to partially radiative modes, the contribution from radiative states is larger than in the 
symmetric waveguide case. 



Figure 11 shows the averaged SE rate F = (2/3)F|| + (1/3)F_l into guided and radiative 



modes for both the symmetric and asymmetric configurations previously analyzed, as a 
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FIG. 9: The spatial dependence of the normaUzed spontaneous emission rate into radiative modes 
for the asymmetric Sihcon waveguide (Si02/Si/air). The same parameters of Fig. |8] have been 
used, (a): the contribution of horizontal dipoles. (b): the contribution of vertical dipoles. 

function of the core thickness. Here Fn = F™ + F™ is the sum over the two polarizations 
for a planar (horizontal) dipole, while F^ = T™ for a vertical dipole and the contributions 
from horizontal (Fy) and vertical (r^) dipoles have been averaged as in the realistic case of 
a randomly oriented dipole in Si. Again, it can be seen that in the asymmetric slab case, 
the contribution of radiative over guided modes in the SE is increased, mainly due to the 
fact that the asymmetric slab supports partially radiative modes that contribute to SE and 
are taken into account explicitly in the present calculation. 

A stronger confinement effect can be achieved in a SOI slot waveguide. The core of 



such a waveguide (see the schematic in Fig. 12) is made up of a very thin layer (slot) of 



low refractive index active material (few tens of nanometers thick) embedded between two 
high-index material regions. In the configuration here considered, the core consists of a 
sequence of Si/Si02 : Er'^+/Si layers and lies on the top of a Si02 cladding grown on a Si 
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FIG. 10: The normalized spontaneous emission rate for the asymmetric Silicon waveguide 
(Si02/Si/air) evaluated as a function of the photonic thickness {ujQd)/c and for a dipole kept 
at the middle of the Silicon core. The same parameters of Figs, [s] and [9] have been used, (a): 
emission from horizontal dipoles. (b): emission from vertical dipoles. For both orientations, the 
onset of emission into guided modes occurs at the threshold value ooqCL/c ~ 0.42. 

substrate. The discontinuity of the normal component of the electric field at the high-index- 
contrast interfaces of the slot gives rise to an increase of the LDoS, which in turn leads to 
an enhancement of SE rate into the waveguide modes. 



In Fig. 13 the calculated SE rates F = (2/3)Fj| + (1/3)F_l into radiative and guided modes for 
a Si Slot waveguide are shown as a function of the dipole position z/\ (even if, in a practical 
case, the Er^"*" emitters are located in the thin Si02 layer). The effect of the discontinuity in 
the z-component of the electromagnetic field at the slot interfaces can be clearly seen: the 
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FIG. 11: The contributions to the total normalized SE rate F = -Tji H — r_|_ from a randomly 

3 3 

oriented dipole into guided and radiative modes, for both symmetric (air/Si/air) and asymmetric 
(Si02/Si/air) Silicon waveguides as a function of the dimensionless thickness lvqcI/c. 



SE rate is mostly due to the decay of vertical (i-oriented) dipoles into TM guided modes [see 
the dashed-dotted line F = F™ in Fig. [T3|^b)], and the total emission into guided modes is 
about six times bigger than the corresponding emission into radiative modes [see the shaded 



regions in Figs. [T3](a) and[T3](b)]. Furthermore, after comparison with Fig. 11, it is evident 
that the light confinement is definitely more effective in a such a Slot waveguide than in a 
symmetric Si waveguide of any core-thickness. 



Moreover, the calculated SE rates shown in Fig. 13 allow to interpret the experimental 
results reported in Ref. 



for the enhancement in the photoluminescence from TM over 
TE polarized modes for a Slot waveguide containing Er'^"'" ions in the oxide (slot) layer. The 



vertical structure is the one depicted in our Fig. [T2| with the same thickness parameters and 
the emission wavelength is 1.54 fim. In the experiment, the TM/TE intensity ratio for light 
emitted from the edge of the waveguide is between 6 and 7.5, with a slight dependence on 



the position of the excitation spot. From Fig. 13, the calculated TM/TE ratio for a dipole 
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Si02 cladding 
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FIG. 12: Schematic of a Slot waveguide. The core consists of a thin slot of Er'^^-doped Si02 having 
thickness ds, embedded between two d2- and d4-thick Si layers; the di-thick Si02 lower cladding is 
grown on a Si substrate. On top of the last Si layer there is air and thus the numerical evaluation 
of the SE rates has been performed for a 6-layers model structure. The following values have been 
used for the layers thicknesses: di=1.9 /.mi, d2=d4=100 nm and d3=20 nm. The values of the 
refractive indices are those which have been used in the structures previously studied. 

embedded in the slot layer is around 7.8 (notice that the r'^^ emission rate is dominated 
by r™, as the TM electric field component in the xy plane has a very small amplitude in 
the slot layer). Thus, the agreement between the theoretical results for the slot waveguide 
obtained within our model and the measurements of Ref. [34j is quite satisfactory. 




VI. CONCLUSIONS 

We have presented a quantum electrodynamical formalism in order to analyze sponta- 
neous emission in generic lossless and non-dispersive multilayer dielectric structures. A 
second quantized form for the electromagnetic field, which follows after its expansion into 
normal modes, has been set up and used to derive the local density of states and express the 
decay rate F as a function of the excited dipole position in the considered structure. The ex- 
pressions derived have been used to study the spontaneous emission in Si-based waveguides 
of different geometries and with realistic parameters. The following conclusions summarize 
our results. 
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FIG. 13: The normalized spontaneous emission rate T = (2/3)r|| + (l/3)r_L for a Slot waveguide 
evaluated as a function of the dipole position, (a): the total emission into radiative modes; the 
refractive index profile is also shown, (b): the total emission into guided modes together with the 
separate contributions of both horizontal dipoles decaying into TE (r™ = T™) and TM modes 
(r™), and vertical dipoles decaying as TM modes only (T™). The layers thicknesses are those 
reported in the caption of Fig. [12] 

(i) The standard basis of radiative states generally used in the description of the elec- 
tromagnetic field modes in a dielectric structure, based on incident /reflected/transmitted 
waves, is not the most appropriate one in radiation emission analysis as it leads to a subtle 
interference between different outgoing components [101 E]. By choosing a set of modes 
specified by a single outgoing radiative component, the total emission rate as well as the 
emission in the upper/lower claddings (more generally, the SE patterns) can be calculated 
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in a simple way, without any interference issue. One basis can be transformed into the other 
after apphcation of the Time-Reversal operation. 

(ii) The evanescent component of partially radiative modes which arise in any asymmetric 
configuration, i.e. when the upper and lower claddings have different refractive indices, is 
relevant for the SE analysis, as it can be seen in the single interface as well as in more 
complicated asymmetric structures. 

(iii) We have calculated and compared SE rates for symmetric (air/Si/air) and asymmetric 
(Si02 /Si/air) silicon waveguides and shown that, in the latter configuration, the lower index 
contrast leads to an increased emission into radiative modes. Such an effect is much more 
evident in a silicon Slot waveguide: in this configuration the discontinuity of the normal 
component of the electromagnetic field which develops at the high-index-contrast interfaces 
of the slot layer, results into an enhancement of the local density of states for TM polarized 
guided modes. As an example, we have analyzed the SE rate in a Si Slot waveguide with 
the same structure parameters used in Ref. |33] and found a very good agreement with 
the experimental evidence of the enhancement of the TM/TE photoluminescence. Thus, 
the model developed turns out to be a useful tool for the analysis of spontaneous emission 
processes in realistic structures such as SOI Slot waveguides. Further work will focus on 
analyzing more complex slot waveguides, as well as photonic crystal slab structures. 
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APPENDIX A: RADIATIVE MODES 



With reference to the geometry of Fig. [T| the field amplitudes of the TE-polarized radia- 
tive modes [see Eqs. ([T]) and ([s]) in Sec. II A| are given by 
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gM+i[XM+ie-^'?M+i(^-^M+i) 



Zj < z < Zj + dj 

Z < Z\ 



(A2) 



(A3) 



For fully radiative modes outgoing in the lower (upper) cladding [see Fig. [2]^b)], Wm+i = 
(Xo = 0) in Eqs. (A1)-(A3) and the amplitude Xq (Wm+i) obtained through the normaliza- 



tion condition Eq. ([3| is given by Xq = (Wm+i = 1/a/^m+i); all the other coefficients 

are then found by application of standard transfer-matrix theory. These results can be 
formally obtained by taking into account a normalization box having width L in the z- 
direction: when L ^ d, d being the thickness of the waveguide core or the thickness of a 
stack of layers in a generic multilayer structure, the contributions from the core/stack are 
of the order 0{d/L) and are negligibly small as compared to the contributions from the 
semi-infinite cladding regions. Thus, the normalization of the radiative modes is determined 
by the cladding regions only, and the values given above are found for the amplitudes Xq and 
Wm+1- When the dielectric constants of the upper and lower claddings are different and the 
conditions for total internal reflection are matched, the modes become partially radiative. 



Without loss of generality, we assume £0 > ^m+i- In this case, when < ku < uj^^^—^, 

c c 

the emission occurs in the lower cladding only and the field becomes evanescent in the upper 
cladding, the z component ^m+i being purely imaginary. The field amplitudes are then found 
through the same conditions given above for the fully radiative modes, together with the 
transformation rule Eq. ([s]), i.e., by taking Wm+i = and replacing gM+i with its complex 



conjugate in Eqs. (Al)-(A3). 



For TM-polarized radiative modes [see Eqs. ^ and (10) in Sec. II A the field amplitudes 
are given by: 
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(A5) 



(A;||,2;) 



QjlZje-'"^'^'-'^-'^^/^^ - y.enA^-^^-'^^l'i)^ , Zj < z < Zj + d^- 



(A6) 



Notice that E™ defined in Eqs. (A5), (A6) have the same dimensions of H"^^ defined in 



Eqs. (A2), (A3), while R™ defined in Eq. (A4) has the same dimensions of E^^"^ defined 



in Eq. (Al). For fully radiative modes outgoing in the lower (upper) cladding, Ym+i = 
(Zo = 0) and the normalization condition Eq. ([s]) yields Zq = 1 (Ym+i = 1) for the 
amplitude of the outgoing component. As for TE-polarized modes, all the other coefficients 
are straightforwardly found after a standard transfer-matrix calculation. For modes which 
are partially radiative in the lower cladding (evanescent in the upper cladding), one takes 
Ym+1 = and replaces ^m+i with its complex conjugate in Eqs. (A4)-(A6). 



APPENDIX B: GUIDED MODES 



The field amplitudes for TE-polarized guided modes [see Eqs. (13) and (14) in Sec. II A 
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(B3) 

where S" is a normalization surface which cancels in the final expressions for the emission 
rates, and the magnetic field is found by application of the Maxwell equation H(r) - 



% c 



Vx 



UJ 



E(r). The M+2 coefficients in the expressions Eqs. (Bl), (B2) and (B3) are found by 



solving the system consisting of M+1 relations which follow from the application of standard 
transfer-matrix theory and the orthormality condition Eq. ([s]) which leads to: 
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with sinc(x) = sin(x)/x. For TM-polarized guided modes [see Eqs. (15) and (16) in Sec. II A 
the field amplitudes are given by 



H™(h,z) 



CM+i/.e-^M+i.^(^-^M+i) ^ z > zm+2 

Dof,e^"^^'-''^ , z < zi 



(B5) 



i 



Z > Zm+1 



h[C,^,e"'^^^'-'^-'^^'^^ + Dj^e-"'^^^'-'^-'^^/^^] , Zj<z< z^ + dj = Zj+i 
Dofye^'-^'-^'K z<zi 



(B6) 



i 



Z > Zm+1 



Bo 



X0m{2-2i) 



Z < Zi 



(B7) 
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where the electric field is obtained from the relation E(r) = — rr^ ^ H(r). As for 

uje{r) 

TE-polarized modes, the M+2 coefficients in the above expressions are derived within the 
transfer-matrix theory together with normalization integral Eq. ([s]) which yields the condi- 
tion: 



H{p,z) p dpdz = 1 



+ 



2xo 2xm 



+1 



\C/ + \D,\')si 



smc 



M 

+ (qD, + D;c,Osinc(fc^ 



{qj - q*) dj 



(B8) 
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